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Abstract. An extended Landau-Lifshitz-Gilbert (LLC) equation is introduced to 
describe the dynamics of inhomogcneous magnetization in a current-carrying wire. 
The coefficients of all the terms in this equation are calculated quantum-mechanically 
for a simple model which includes impurity scattering. This is done by comparing 
the energies and lifetimes of a spin wave calculated from the LLC equation and from 
the explicit model. Two terms are of particular importance since they describe non- 
adiabatic spin-transfer torque and damping processes which do not rely on spin-orbit 
coupling. It is shown that these terms may have a significant influence on the velocity 
of a current-driven domain wall and they become dominant in the case of a narrow 
wall. 
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1. Introduction 

The effect of passing an electric current down a ferromagnetic wire is of great current 
interest. If tfie magnetization is inliomogeneous it experiences a spin-transfer torque 
due to tlie current [H El El H]. The effect is described phenomenologically by adding 
terms to the standard LLG equation [5l E]. The leading term in the spin-transfer 
torque is an adiabatic one arising from that component of the spin polarization of the 
current which is in the direction of the local magnetization. However , in considering 
the current-induced motion of a domain wall, Li and Zhang [3l H] found that below a 
very large critical current the adiabatic term only deforms the wall and does not lead 
to continuous motion. To achieve this effect they introduced [7] a phenomenological 
non-adiabatic term associated with the same spin non- conserving processes responsible 
for Gilbert damping. Subsequently Kohno et al [8] derived a torque of the Zhang-Li 
form quantum-mechanically using a model of spin-dependent scattering from impurities. 
This may arise from spin-orbit coupling on the impurities. More recently Wessely et 
al [9] introduced two further non-adiabatic terms in the LLG equation in order to 
describe their numerical calculations of spin-transfer torques in a domain wall. These 
quantum-mechanical calculations using the Keldysh formalism were made in the ballistic 
limit without impurities and with spin conserved. Other terms in the LLG equation, 
involving mixed space and time derivatives, have been considered by Sobolev et al |12j . 
Tserkovnyak et al ^lOj, Skadsen et al [11] and Thorwart and Egger |13j . 

The object of this paper is to give a unified treatment of all these terms in the LLG 
equation and to obtain explicit expressions for their coefficients by quantum-mechanical 
calculations for a simple one-band model with and without impurity scattering. The 
strategy adopted is to consider a uniformly magnetized wire and to calculate the effect 
of a current on the energy and lifetime of a long wavelength spin wave propagating along 
the wire. It is shown in section 2 that coefficients of spin-transfer torque terms in the 
LLG equation are directly related to q and q^ terms in the energy and inverse lifetime 
of a spin wave of wave- vector q. The Gilbert damping parameter is the coefficient of the 
uj term in the inverse lifetime, where uj is the spin-wave frequency. It corresponds to the 
damping of a g = spin wave while higher order terms ujq and ujq^ relate to damping 
of spin waves with finite wave-vector q. The relation between the q term in the spin 
wave energy and the adiabatic spin-transfer torque has been noticed previously j2l [T^. 
We find that the q term in the spin wave lifetime relates to the Zhang-Li non-adiabatic 
spin transfer torque. Our result for the coefficient of the Zhang-Li term is essentially the 
same as that obtained by Kohno et al [8] and Duine et al [15] but our derivation appears 
simpler. The q^ terms in the spin wave energy and lifetime are related to the additional 
non-adiabatic torques we introduced into the LLG equation [9], together with an extra 
one arising from spin non-conserving scattering. Explicit expressions for the coefficients 
of these terms are obtained in section 3. In section 4 we discuss briefly the importance 
of the additional terms in our extended LLG equation for current-driven motion of a 
domain wall. Some conclusions are summarized in section 5. 
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2. The LLG equation and spin waves 



We write our extended LLG equation in the dimensionless form 
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Here s(z, t) is a unit vector in the direction of the local spin polarisation, time t is 
measured in units of (7/io^s)~^ and the coordinate z along the wire is in units of the 
exchange length l^x = (2A/yUo?Ti^)^/^. The quantities appearing here are the gyroscopic 
ratio 7 = 2/iB/^,the permeability of free space /xq and two properties of the ferromagnetic 
material, namely the saturation magnetisation and the exchange stiffness constant 
A. is a unit vector in the z direction along the wire. The equation expresses the 
rate of change of spin angular momentum as the sum of various torque terms, of which 
the «!, a'l, a, /, ai, /i and gi terms are proportional to the electric current flowing. 
The second term in the equation is the standard Gilbert term, with damping factor 
a, while the a[ and ctg terms introduce corrections for spin fluctuations of finite wave- 
vector. Skadsem et al [H] point out the existence of the a'2 term but do not consider 
it further. It was earlier introduced by Sobolev et al |12] within a microscopic context 
based on the Heisenberg model. The ai and a2 terms are found to renormalise the spin 
wave frequency, but for the model considered in section 3 we find that Oi is identically 
zero. We shall argue that this result is model- independent. Tserkovnyak et al [10] and 
Thorwart and Egger [13] find non-zero values of ai which differ from each other by a 
factor 2; they attribute this to their use of Stoner-like and s — d models, respectively. 
Thorwart and Egger [13] also find the a'-^ term and they investigate the effect of ai and a'l 
terms on domain wall motion. Their results are difficult to assess because the constant 
|s| = 1 is not maintained during the motion. In eq.([l]) we have omitted terms involving 
the second order time derivatives, whose existence was pointed out by Thorwart and 
Egger [13]; one of these is discussed briefly in section 3.2. 

The first term on the right-hand side of eq. ([T]) is due to exchange stiffness and the 
next term arises from an external magnetic field B^xt^z with dimensionless coefficient 
bext = Bext/ fJ'O^s- The third term is the adiabatic spin transfer torque whose coefficient 
a is simple and well-known. In fact [Sj H] 
1 hJP 

a = ^ ^ (2) 
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where J is the charge current density and e is the electron charge (a negative quantity). 
The spin polarisation factor P = (J| — J\)/{J^ + Jj), where J|, J| are the current 
densities for majority and minority spin in the ferro magnet (J = + J|). Eq.(l2]) is 
valid for both ballistic and diffusive conduction . The fourth term on the right-hand 
side of eq.dl]) is the Zhang-Li torque which is often characterised [H] by a parameter 
(3 = f /a. The next term is the Ei term of eq.(7) in ref. [9]. It is a non-adiabatic torque 
which is coplanar with s{z) if s{z) lies everywhere in a plane. As shown in ref. [9] it 
is the z derivative of a spin current , which is characteristic of a torque occurring from 
spin-conserving processes. In fact this term takes the form 
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The /i term may be written in the form 
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If s{z) lies in a plane, the case considered in ref. [9J, the first term vanishes and we 
recover the Fi term of eq.(9) in ref. [9]. Its derivative form indicate that it arises from 
spin-conserving processes so we conclude that the coefficient /i is of that origin. This 
is not true of the last term in eq.([T]) and we associate the coefficient gi with spin non- 
conserving processes. For a spin wave solution of the LLG equation, where we work 
only to first order in deviations from a state of uniform magnetisation, the last three 
terms of eq.([T]) may be replaced by the simpler ones 

s s 

Apart from additional terms, eq.([T]) looks slightly different from eq.(7) of ref. ^ because 
we use the spin polarisation unit vector s rather than the magnetisation vector m and 
s = — m. Furthermore the dimensionless coefficients will take different numerical values 
because we have used different dimensionless variables z and t to avoid introducing the 
domain wall width which was specific to ref. [9]. The torques due to anisotropy fields 
were also specific to the domain wall problem and have been omitted in eq.([T]). 

We suppose that the wire is magnetised uniformly in the z direction and consider a 
spin wave as a small transverse oscillation of the spin polarisation about the equilibrium 
state or, when a current fiows, the steady state. Thus we look for a solution of eq.([I]) 
of the form 

s = (ce^^""-"*), de'^'''-^'\ -1) (6) 

where the coefficients of the x and y components satisfy c <^ 1, c? ^ 1. This represents 
a spin wave of wave-vector q and angular frequency uo propagating along the z axis. 
When (El) is substituted into eq.([I]) the transverse components yield, to first order in c 
and d, the equations 

— iXc + fid = , /ic + iXd = (7) 
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where 

\ = u — aq + aiq^ — a2Ujq^ + ia'^qu 

fi= - iau + bext + + ifq + i{fi + gi)q^ + aiujq - iuq^a^- (8) 

On eliminating c and d from eq.(l7]) we obtain = /i^. To obtain a positive real part 
for the spin wave frequency, we take A = /x. Hence 

uj[l- aiq - a2q^) 
= bext + aq + q^ - aiq^ 

+ % {uj (-« - a'lg - ^ + ^ ^ _ 

Thus the spin wave frequency is given by 

ijj = iuj2 (10) 

where 

cui ~ (l - aiq - aag^) ^ [b^xt + aq + q^ - aiq^) 

u;2 - (l - otiq - ot2q^) ^ [t^i (a + + «29^) - fQ - (/i + fl'i) 9^] • (H) 

Here we have neglected terms of second order in a, a[, a'2, f, fi and gi, the coefficients 
which appear in the spin wave damping. This form for the real and imaginary parts 
of the spin wave frequency is convenient for comparing with the quantum-mechanical 
results of the next section. In this way we shall obtain explicit expressions for all the 
coefficients in the phenomenological LLG equation. Coefficients of odd powers of q are 
proportional to the current flowing whereas terms in even powers of q are present in the 
equilibrium state with zero current. 



3. Spin wave energy and lifetimes in a simple model 

As a simple model of an itinerant electron ferromagnet we consider the one-band 
Hubbard model 

Ho = -t^^ c\„Cja + rii^^riii - fisBext (^it - n^) , (12) 

ija i i 

where creates an electron on site i with spin a and = cl^Cia- We consider a simple 
cubic lattice and the intersite hopping described by the first term is restricted to nearest 
neighbours. The second term describes an on-site interaction between electrons with 
effective interaction parameter U; the last term is due to an external magnetic field. It 
is convenient to introduce a Bloch representation, with 

^ic- = ;;^5Z^''''"^l ' ^k. = 4ck., (13) 

* i 

Ck = —t e^^ ''^ = — 2t (cos kxttQ + cos kya^ + cos fczOo) . (14) 
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The sum in eq.( lT^ is over all lattice cites Rj whereas in eq. (fT^ pi = 
(±ao, 0, 0), (0, ±ao, 0), (0, 0, ±ao) are the nearest neighbour lattice sites. Then 

Hq = ^ Ck^ka + rii^riii - fisBext ^ (^kt - ^ki) • (15) 

kcr j k 

To discuss scattering of spin waves by dilute impurities we assume that the effect of 
the scattering from different impurity sites adds incoherently; hence we may consider 
initially a single scattering center at the origin, We therefore introduce at this site a 
perturbing potential u + vl ■ a, where 1 = (sin^ cos0, sin(/) sin ^, cos^) is a unit vector 
whose direction will finally be averaged over, u is the part of the impurity potential 
which is indepndent of the spin cr and the spin dependent potential f 1 ■ a is intended to 
simulate a spin-orbit L-a interaction on the impurity. It breaks spin rotational symmetry 
in the simplest possible way. Clearly spin-orbit coupling can only be treated correctly 
for a degenerate band such as a cZ-band, where on-site orbital angular momentum L 
occurs naturally. The present model is equivalent to that used by Kohno et al [8] and 
Duine et al p^. In Bloch representation the impurity potential becomes V = Vi + V2 
with 

kik2 kik2 

V2 = ve-^'^sme^ it^k^i + ve'^sme^ ^1 ii^^^T (16) 

kik2 kik2 

and vi = u + vcos9 , v^^ = u — vcosO. To avoid confusion we note that the spin 
dependence of the impurity potential which occurs in the many-body Hamiltonian Hq+V 
is not due to exchange, as would arise in an approximate self consistent field treatment 
(e.g. Hartree-Fock) of the interaction ?7 in a ferromagnet. 

3.1. Spin wave energy and wave function 

In this section we neglect the perturbation due to impurities and determine expressions 
for the energy and wave function of a long-wave length spin wave in the presence of an 
electric current. The presence of impurities is recognised implicitly since the electric 
current is characterised by a perturbed one-electron distribution function /ka which 
might be obtained by solving a Boltzman equation with a collision term. We consider 
a spin wave of wave-vector q propagating along the z axis, which is the direction 
of current flow. Lengths and times used in this section and the next, except when 
specified, correspond to actual physical quantities, unlike the dimensionless variables 
used in section 2. 

We first consider the spin wave with zero electric current and treat it, within 
the random phase approximation (RPA), as an excitation from the Hartree-Fock (HF) 
ground state of the Hamiltonian f|T5|) . The HF one electron energies are given by 

Eka = ek + U {n-a) - fJ.B<TBext (17) 
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where a = 1,-1 for | and | respectively, and {n_„) is the number of —a spin electrons 
per site. In a self-consistent ferromagnetic state at T = 0, {n^) = /ko- and 

n = ^^{ua), where, /ko- = 0{Ef — E^^), n is the number of electrons per atom, and 
Ep is the Fermi energy. is the number of lattice sites and 9{E) is the unit step 
function. The spin bands i^kcr given by eq. (|T711 are shifted relative to each other by an 
energy A + where A = U{n^ — n^) is the exchange splitting. The ground state 

is given by |0) = Hkcr'^ko-I) "where |) is the vacuum state and the product extends over 
all states kcr such that /ko- = 1- Within the RPA, the wave function for a spin wave of 
wave-vector q, excited from the HF ground state, takes the form 

|q) = ^qE^kcUqi^kT|0) (18) 

k 

where A'^q is a normalisation factor. The energy of this state may be written 

Ed = Egr + huJ^ = Egr + 2flBBext + ^^^q (19) 

where Egr is the energy of the HF ground state and hw^ is the spin wave excitation 
energy. On substituting ( |T8i) in the Schrodinger equation {Hq — -Eq)|q) = and 
multiplying on the left by (0|cj^,-|.Ck'-q|, we find 

Ak' (ek'+q -e^, + A-fhj'^)=^Yl ^k/kT (1 - /k+qi) . (20) 

k 

Hence we may take 

Ak = A (ek+q -e^ + A-hw'^) ~' (21) 
and, for small q, hu'^ satisfies the equation 

I ^ ^ /kt - /k+qj ^22) 

iV^ek+q-ek + A-;Lu;' 

This is the equation for the poles of the well-known RPA dynamical susceptibility x(q, ^) 
[TB] . The spin wave pole is the one for which hw'^ — as q ^ 0. 

To generalise the above considerations to a current- carrying state we proceed as 
follows. We re-interpret the state |0) such that (0| . . . |0) corresponds to a suitable 
ensemble average with a modified one-electron distribution /ko-- When a current fiows 
in the z direction we may consider the | and | spin Fermi surfaces as shifted by small 
displacement ^^k^, ^^k^ where k^ is a unit vector in the z direction. Thus 

/ko = 0{Ef - ^k+5,L,o) 

- e{EF - Eko) - SJ{Ef - ^ko)^ (23) 
and the charge current density carried by spin a electrons is 
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where {{de]^/dkzY)„ is an average over the a spin Fermi surface and p^{Ep) is the 
density of a spin states per atom at the Fermi energy, 
following related quantities; 



We shall also encounter the 




(25) 



(26) 

To derive eqs. (12^ and 5^ has been eliminated using eq. (121l) . 

To solve eqn. fl22|) for hjj'^ we expand the right-hand side of the equation in powers 
of (^k+q — Ck — ^q) / ^ ^"^^ make the further expansions 
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(27) 

Bq + Dq' + Eq' + ... (28) 

in powers of q. We retain all terms up to q^ except those involving B^\ the coefficients 
B and E are proportional to the current and we keep only terms linear in the current. 
Hence we find a solution of eq. (l22|) in the form fl28l) with 
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Here is the total number of a spin electrons so that A^^. = N{n^). 

In the absence of spin-orbit coupling the expression for B in terms of spin current 
is a general exact result even in the presence of disorder, as shown in Appendix A. The 
coefficient D is the standard RPA spin- wave stiffness constant (e.g ref. [16]). We note 
that, in the limit A cxd, E takes the simple form —aQB/6. 

On restoring the correct dimensions (as indicated after eq.([T])) to the expression 
for Ui in eq.l fTTl) we may determine the coefficients a and ai by comparing with the 
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equation 



hu^ = 2^iBBext + Bq + Dq^ + Eq^ 



(32) 



From the coefficient of q we have 

a + aibext = B/ (2/iB/io"^Jez) • 



(33) 



a and B are both determined directly from the spin current JP independently of a 
particular model (see appendix A) so that bext should not enter their relationship. We 
conclude quite generally that ai = 0. In this case we find that on combining eqs. (l33!) 
and (!29ll . and noting that = —fiB{Ni — Ni)/Nal, eq.(l2]) is obtained as expected. 
In section 3.2 we show explicitly for the present model that ai = 0. This conflicts 
with the results of refs. [TO] and [13]. From the coefficients of in eqs. ffTT]) and fl32|) we 
find 1 + a2bext = D / {A^bA/itls). Thus an external field slightly disturbs the standard 
relation A = DitIs/A^b- However in the spirit of the LLG equation we take A and 
which enter the units of length and time used in eq.([l]), to be constants of the 
ferromagnetic material in zero external field. The coefficients of q^ in eqs. fllll) and fl32l) 
yield the relation (taking ai = 0), 



We defer calculation of 02 until section 3.2 and the result is given in eq.( H^ . Combining 
this with eqs. (lMI) and (!3T!) we find 



We have thus derived an explicit expression , for a simple model, for the coefficient ai 
of a non-adiabatic spin torque term which appears in the LLG equation ([T]). We have 
neglected the effect of disorder due to impurities . In the absence of spin-orbit coupling 
the expression for the adiabatic torque coefficient a, given by eq.([2]), is exact even in 
presence of impurities. In the next section we shall calculated further non-adiabatic 
torque terms, with coefficients /i and gi, as well as damping coefficients a, a'^ and a^- 
In the present model all these depend on impurity scattering for their existence. 

3.2. Spin wave lifetime 

The solutions of eq. (l22l) are shown schematically in figure [H They include the spin 
wave dispersion curve and the continuum of Stoner excitations c|^_,_q|Ckt |0) with energies 
-f'k+qi — -E'kt • The Zeeman gap 2fiBBext in the spin wave energy at q = does not appear 
because we have plotted hu'^ rather than (see eq. (fT9!) ). Within the present RPA the 
spin wave in a pure metal has infinite lifetime outside the continuum and cannot decay 
into Stoner excitations owing to conservation of the momentum q. However, when the 
perturbation Vi due to impurities is introduced (see eqn. f|T6l) ). crystal momentum is no 
longer conserved and such decay processes can occur. These are shown schematically 
by the dotted arrow in figure [H If the bottom of the [ spin band lies above the Fermi 
level there is a gap in the Stoner spectrum and for a low energy (small q) spin wave 



ai + = E/ (2fiBfJ'0'msllx) ■ 



(34) 




(35) 
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such processes cannot occur. However the spin-flip potential V2 enables the spin wave 
to decay into single particle excitations cl.^^^Cka\0) about each Fermi surface and these 
do not have an energy gap. 

The lifetime of a spin wave of wave- vector q is thus given simply by the "golden 
rule" in the form 

27T 



Tq ^ = -^Nimp (Ti + T2) 



(36) 



where iVj„„ is the number of impurity sites and 
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(37) 



We first consider Ti and, using eqns. (fT6!) and (ITS!) . we find 
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q) = ^/kT (1 - fpl) l^kVi (1 - fpl) - Ap-qt^T/p-q?] 
= ^/kT(l-/pi)(^kn-^p-q^T) 



(3^ 



for small q. The last line follows from two considerations. Firstly, because of the S- 
function in eg. (1371) we can consider the states and P| to be close to their respective 
Fermi surfaces. Secondly the | spin Fermi surface lies within the t Fermi surface and q 
is small. Hence 

^1 = ]^ E (1 - /pi) ^ (^q - Epl + ^kt) (^k^ - A^^^v^f . (39) 



kp 
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To evaluate this expression in the case when a current flows we use the distribution 
function /k^- given by eq. fl23l) . Thus, neglecting a term proportional to the square of the 
current, we have 

kp 



X 



e{EF - ^kT)^(^pi - Ep) - 6^e{Ep^ - Ef)5{Ef - E^{ 



de ' 



(40) 



We wish to expand this expression, and a similar one for T2, in powers of q to 0{q^) so 
that we can compare with the phenomeno logical expression (eq.l fTTl) ) for the imaginary 
part of the spin wave frequency, which is given by r~^/2. It is straight-forward to expand 
the second factor in the above sum by using eqs. fl2T]) and fl28|) . We shall show that the 
contribution to Ti of the flrst term in square brackets in eq. fHUl) leads to a contribution 
proportional to spin wave frequency tUq. Together with a similar contribution to T2 it 
yields the Gilbert damping factor a as well as the coefficients a[, a'^ of the terms in 
eq. ffTT]) which give the q dependence of the damping. The remaining terms in eq. fHOl) 
yield the spin-transfer torque coefficients /, j\ and g\. 

The normalisation factor which appear in eq. lHOl) leads naturally to the factor 
(1 — a\q — a2q'^)~^ which appears in eq. ffTT]) . From eq. lITS!) it is given by 

§ E (^^^T - K-M ■ (41) 
in powers of q, and using eq. fl2^ . we find to O(g^) that 

[e{EF - ^kt) - 0{Ef - Eki)] } . (42) 

^- {ly-i - lyi) — \uf.^/ 
We deduce that 

tti = (43) 




and 

a2 



-) E (||) ' [^(^^ - ^"^t) - 0{Ef - i?ki)] . (44) 



The result ai = 0, which was predicted on general grounds in section 3.1 and in 
Appendix 1, arises here through the absence of a g term, proportional to current, 
in the spin wave normalisation factor. In the derivation of eq. fH2|) this occurs due 
to a cancellation involving the Bq terms in the spin energy, which appears in Ak. 
Without this cancellation we would have ai = 2B/l^r^/S. which is of the form obtained 
by Tserkovnyak et al [10] and Thorwart and Egger [T3] . 
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We now return to the programme for calculating the LLG coefficients 
a, a[, a'2, f, fi, Qi which was outlined after eq. fHOj) . We have seen that the q dependence 
of A*"^ corresponds to the prefactor in eq.f fTTl) . Hence to determine the coefficients listed 



above we can take = Nq = {N-^ — N^) in Ti and T2 when we expand terms in 



powers of q to substitute in eq.( ]36l) and compare with eq. ljTTll . We ffist consider the 
case g = in order to determine the Gilbert damping factor a. Thus only the ffist term 
in square brackets in eq. fHOl) contributes, since de^/dkz is an odd function k^, and 



Ti(q = 0) 



4f ^cos^ 9 



kp 



^kt) 0{Ep - E^^)d{Ep^ - Ef) (45) 



where cos^ 6 is an average over the angle appearing in the impurity potential V (eq. f|T6l) ) 

0. The summations in eq. 



and we shall assume cos 6' = 0. The summations in eq. fBS]) may be replaced by energy 
integrals involving the density of states of per atom Po-(e) of the states i^ko-- Then, to 
order (^q)^; 



Ti(q = 0) 



4f ^cos^ 6 



where Pcr(e) and its derivative p^(e) are evaluated at e 



Ep. Similarly 



T2(q = 0) 



f ^sin^ 9 



Nf-Ni 



hujQ (p^ + pf) 



(46) 



(47) 



and no terms appear. We have included the uJq term in eq. fl46l) merely because it 
corresponds to a term s x x |^ j in the LLG equation whose existence was noted by 
Thorwald and Egger [13] . We shall not pursue terms with second-order time derivatives 
any further. Since the imaginary part of the spin wave frequency is given by t~^/2 it 
follows from eqs.dH]), ([36]), (j46]) and (j!?!) that 



TTCV 



a 



4cos2 Op-^pi + sin^ 6 (p^ + pf) 



(48) 



where c = Nimp/N is the concentration of impurities, in agreement with Khono et al [8] 
and Duine et al [1^. If the direction of the spin quantisation axis of the impurities is 
distributed randomly cos^ 6 = 1/3, sin^ 6* = 2/3 so that a is proportional to (pt + Pi)^- 
To investigate the q dependence of Gilbert damping, and thus evaluate a'l and a'2 
in eq. fllll) . the second factor in the summation of eq. fllOl) must be expanded in powers 
of q. All the terms which contribute to the sum are of separable form g{]s.)h{p). The 
contribution to Ti of interest here , proportional to cUq, again arises from the ffist term 
in square brackets in eq. (l40]) . and similarly for T2. The summations required in eq. (l40!) 
are of the form 



Y,H^n~ ^Pi + ^kt) 0{Ef - E^^)d{Ep^ - EF)g{\^)h{p) 
kp 



(^(k))t (/i(k))| pipi^^, 



(49) 
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where ((?(k))o- = ^ J2]i90^)^i^F — -^ko-) is an average over the Fermi surface, as used 
previously in section 3.1. After some algebra we find 

a[ = 25a/AZex (50) 



TTC 



an 




+ 



f ^sin^ 9 
2Da 



dk'i 




Ml 



(51) 



We note that, unlike a and a'^, the coefficient a'2 is non-zero even when the spin- 
dependent part of the impurity potential, f , is zero. In this case the damping of a 
spin wave of frequency uj and small wave-vector q is proportional to p^piu^uq^. In 
zero external field a; ~ so that the damping is of order g^. This damping due to 
spin-independent potential scattering by impurities was analysed in detail by Yamada 
and Shimizu [TTj. One of the Fermi surface averages in eq. (l5Tl) is easily evaluated using 
eqs.dllD and Thus 



(52) 



In the spirit of the LLG equation we should take Bext = in evaluating the coefficients 

We now turn to the evaluation of the non-adiabatic spin-transfer torque coefficients 
/, /i and gi. These arise from the second and third terms in square brackets in eq.l HOl) . 
and in a similar expression for T2. The summations involved in these terms differ from 
those in eq. (H9il since one 6'-function is replaced by a 5-function. This leads to the 
omission of the frequency factor The Fermi surface shifts are elininated in 

favour of currents Jo- by using eq. fl2^ . 

By comparing the coefficient of q in the expansion of eq. (l36i) with that in eq. lfTTil 
we find the coefficient of the Zhang-Li torque in the form 

•2 h 



f 



TTCV 



2cos2 6 (pt Jj, - Pi J|) + sin 6 (p| J| - p^J^] 



Pom^Alex e 

This is in agreement with Khono et al [S] and Duine et al 



(53) 



In the "isotropic" 



impurity case, with cos^ 6 

a U{p^ + pi) 



P 



1/3, sin^ e = 2/3, it follows from eqs.([53]), (gH]) and that 
2 



(54) 



In the limit of a very weak itinerant forromagnet p„ p, the paramagnetic density of 
states, and Up — 1 by the Stoner criterion. Thus in this limit (3 = a. Tserkovnyak et 
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al [TU] reached a similar conclusion. For a parabolic band it is straightforward to show 
from Stoner theory that j3/a > 1 and may be as large as 1.5. 

As discussed in section 2 the coefficient /i is associated with spin- conserving 
processes, and hence involves the spin independent potential u. The coefficient gi 
is associated with spin non-conserving processes and involves v. By comparing the 
coefficient of in the expansion of eq. (|36|) with that in eq. ffTTl) we deduce that 



and 



/i 



91 



TTC 



2/^0^2/3 



V {Ki + 2Li + Ml] 



(55) 



''ex 
TTCf 2 



3D 



6 



/ 



2/iom2/3^ 



cos2 e {5Ki + 6Li - Ml) + sin' ^ (3^2 + 4L2) 



(56) 



Here 



Ki 

Ml 




^2 = Kipi + K^p^ 
L2 = Lipi — L^p^ 



+ A 



5V 



(57) 



This complete the derivation of expressions for all the LLG coefficients of eq.([T]) within 
the present impurity model 



4. The extended LLG equation applied to current-driven domain wall 
motion 

In a previous paper [9j we introduced the ai and /i terms of the extended LLG equation 
(cf. eqns.([T]), ([3]) and (j4])) in order to describe numerically-calculated spin-transfer 
torques acting on a domain wall when it is traversed by an electric current. In that 
work the origin of the small fi term for a pure ferromagnetic metal was specific to 
the domain wall problem; it was shown to be associated with those electronic states at 
the bulk Fermi surface which decay exponentially as they enter the wall. The analytic 
derivation of /i in section 3 (see eqn.( |55ll ) is based on impurity scattering in the bulk 
ferromagnet and applies generally to any slowly-varying magnetization configuration. 
For a ferromagnetic alloy such as permalloy both mechanisms should contribute in the 
domain wall situation but the impurity contribution would be expected to dominate. 

To describe a domain wall we must add to the right-hand side of eqn.([T]) anisotropy 
terms of the form 

— (s ■ e^) s X ej, + b^^ (s ■ e^) s x e^, (58) 

where e^ is a unit vector perpendicular to the plane of the wire. The first term 
corresponds to easy-plane shape anisotropy for a wire whose width is large compared 
with its thickness and the second term arises from a uniaxial field Hu along the wire, 
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so that b = nis/Hu- The solution of eqn.([T]), with the additional terms for a 

stationary Neel wall in the plane of the wire, with zero external field and zero current, 
is 

s = (sec/i(z/6^/2), 0, -tanh{z/b^/^)). (59) 

As pointed out in ref. [9J there is no solution of the LLG equation of the form 
s = F{z — vwt), corresponding to a uniformly moving domain wall, when the /i term 
is included. It is likely that the wall velocity oscillates about an average value, as 
predicted by Tatara and Kohno fiSl [T9] for purely adiabatic torque above the critical 
current density for domain wall motion. However, we may estimate the average velocity 
vw using the method of ref. [9j. The procedure is to substitute the approximate form 
s = F{z — Vwt) in the extended LLG equation ([T]), with the terms (l58l) added, take the 
scalar product with F x F' and integrate with respect to z over the range (— oo, oo). 
The boundary conditions appropriate to the wall are s as z ±oo. Hence for 

bext = we find the dimensionless wall velocity to be 

_ / IZ (F X Ff dz + h JZ (F X F-)^ dz + g, JZ (F'f dz 

aJZiF^Ffdz + a'^JZiF'fdz ' ^'"^ 

To estimate the integrals we take F{z) to have the form of the stationary wall s{z) 
(eqn. (l59!l ) and, with the physical dimensions of velocity restored, the wall velocity is 
given approximately by 

(3 1 + f,{3fb)-' 
a 1 + a2{ab) ^ 

where vq = hbPJ I ijnse). We have neglected gi here because, like / and a, it depends 
on spin-orbit coupling but is a factor (do/^ex)^ smaller than / (cf. eqns.( l53|) and (!56|) ). 
/i and a'2 are important because they do not depend on spin-orbit coupling. 

It is interesting to compare vw with the wall velocity observed in permalloy 
nanowires by Hayashi et al |20j. We first note that vq is the velocity which one obtains 
very simply from spin angular momentum conservation if the current-driven wall moves 
uniformly without any distortion such as tilting out of the easy plane and contraction 
|21] . This is never the case, even if /i = 0, = 0, unless (3 = a. For a permalloy 
nanowire, with ^^ms = 1 T, i;o = HOP m/s for J = L5 ■ lO^A/cm^. Thus, from the 
standard theory with /i = 0, = 0, vw = llOPP/a m/s for this current density. In 
fact Hayashi et al [20] measure a velocity of 110 m/s which implies (3 > a since the spin 
polarization P is certainly less than 1. They suggest that f3 cannot exceed a and that 
some additional mechanism other than spin-transfer torque is operating. However in the 
discussion following eqn. flS^ we pointed out that in the model calculations it is possible 
to have (5 > a. Even if this is not the case in permalloy we can still have vw > vo if the 
last factor in eqn. (l6T]) is greater than 1 when /i and are non-zero. We can estimate 
terms in this factor using the observation from ref. [20] , that Iw = hxb^^^ = 23 nm, where 
Iw is the width of the wall. From eqns. (!53i) and (!55ll we find /i/(/6) ~ {u/v)^{kFlw)~^, 
where kp is a. Fermi wave- vector. In permalloy we have Fe impurities in Ni so that in 
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the impurity potential u + v -a we estimate m ~ 1 eV and v ~ 0.005 eV. The value for v 
is estimated by noting that the potential v ■ a is intended to model spin-orbit coupling 
of the form • a with ^ < 0.1 eV and (L^) pe ~ 0.05, Lz being the component of orbital 
angular momentum in the direction of the magnetization [22]. Hence w/f ~ 200 and 
kplw ~ 200 so that fi/{fb) ~ 1. ay^ab) is expected to be of similar magnitude. We 
conclude that the a2 and /i terms in the LLG equation ([1]) can be important in domain 
wall motion and should be included in micromagnetic simulations such as OOMMF 
[23] . For narrower domain walls these terms may be larger than the Gilbert damping a 
and non-adiabatic spin- transfer torque / terms which are routinely included. Reliable 
estimates of their coefficients are urgently required using realistic multiband models of 
the ferromagnetic metal or alloy. 

5. Conclusions 

The coefficients of all the terms in an extended LLG equation for a current-carrying 
ferromagnetic wire have been calculated for a simple model. Two of these (/i and 
are of particular interest since they do not rely on spin-orbit coupling and may 
sometimes dominate the usual damping and non-adiabatic spin-transfer torque terms. 
One term (ai) which has been introduced by previous authors is shown rigorously to 
be zero, independent of any particular model. Solutions of the extended LLG equation 
for domain wall motion have not yet been found but the average velocity of the wall is 
estimated. It is pointed out that the fi and ^2 terms are very important for narrow 
walls and should be included in micrmagnetic simulations such as OOMMF. It is shown 
that there is no theoretical reason why the wall velocity should not exceed the simplest 
spin-transfer estimate vq, as is found to be the case in experiments on permalloy by 
Hayashi et al [20] 
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Appendix A. 

The simple single-band impurity model used in the main text is useful for obtaining 
explicit expressions for all the coefficients in the LLG equation ([1]). Here we wish to 
show that some of these results are valid for a completely general system. We suppose 
the ferromagnetic material is described by the many-body Hamiltonian 



H = Hi + Hint + Hext 



(A.l) 



where Hi is a one-electron Hamiltonian of the form 



Hi=Hk + Hso + V. 



(A.2) 
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Here Hj^ is the total electron kinetic energy, Hgo is the spin-orbit interaction, V is a. 
potential term. Hint is the coulomb interaction between electrons and H^xt is due to an 
external magnetic field Bext in the z direction. Thus 

Hext = —'^fJ^sSQEext 

where 5*° is the z component of total spin. Both Hgo and V can contain disorder. Since 
we are interested in the energy and lifetime of a long-wavelength spin wave we consider 
the spin wave pole, for small q, of the dynamical susceptibility. 



X(q,c.) = / dt{{S^it),S^^))e-'-'' (A.4) 

density. Here 



[u- = uj — ie) where = iS^ are Fourier components of the total transverse spin 



{{S-{t),S^^)) = ^{[S-{t),Sl^])9it). (A.5) 

In general we shall take the average () in a steady state in which a charge current density 
J is fiowing in the q direction. Following the general method of Edwards and Fisher 
we use equations of motion to find that 

where hbext = 2i^BBext, = [S^,Hi] and 

Xciq^oo) = J dt{{C^{t),C\))e-'^\ (A.7) 
For small q and tu, x is dominated by the spin wave pole, so that 

x{<iM = -^r-^ — T (^-^^ 

where bext + is the spin wave frequency, in general complex corresponding to a finite 
lifetime. Following ref. p4J we compare (1A.6P and flA.SP in the limit ojq ^ oj — bext to 
obtain the general result 



Edwards and Fisher [24] were concerned with Re Ug whereas Kambersky [25] derived the 
above expression for ImuOq for the case g = 0, and zero current fiow. His interest was 
Gilbert damping in ferromagnetic resonance. Essentially the same result was obtained 
earlier in connection with electron spin resonance, by Mori and Kawasaki [2S], see also 
Oshikawa and Affleck [27]. Since 5*" commutes with the potential term V, even in the 
presence of disorder, we have 

= [S~, H{\ = [S-, Hk] + [S-, Hso] . (A. 10) 

For simplicity we now neglect spin-orbit coupling so that 

= [S^,Hk\ =hqJ- (A.ll) 
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where the last equation defines the spin current operator . For a general system, with 
the n^^ electron at position r„ with spin and momentum p„, 



s: 



Hence, from eqns.( IX?TTl) and ( 1A.12I) . 



2m 



(A.12) 



(A.13) 



where is the total number of electrons and = Pn/i^ is the electron velocity, so 
that eJ2n('^n^n) is the total spin current. Hence from eq. (lA.9p . we find 



with 



B 



hm Xj{0,^) 



2{Sq) [2m i^-'bex 



Bq 



(A. 14) 



hfiBPJ/ems. (A.15) 

This expression for B has been obtained by Bazaliy et al [2] and Fernandez-Rossier et 
al [H] for simple parabolic band, s — d and Hubbard models. The derivation here is 
completely general for any ferromagnet, even in the presence of disorder due to impurities 
or defects, as long as spin-orbit coupling is neglected. Eqs. ([2]) and flA.MD are both valid 
for arbitrary fogxt, so that in eq. ( !33l) we must have ai = 0. 
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